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We study the scaling limit of a fully packed loop model in two dimensions, where the loops are 
endowed with a bending rigidity. The scaling limit is described by a three-parameter family of 
conformal field theories, which we characterize via its Coulomb-gas representation. One choice for 
two of the three parameters reproduces the critical line of the exactly solvable six-vertex model, 
. while another corresponds to the Flory model of polymer melting. Exact central charge and critical 

exponents are calculated for polymer melting in two dimensions. Contrary to predictions from 
mean-field theory we show that polymer melting, as described by the Flory model, is continuous. 
We test our field theoretical results against numerical transfer matrix calculations. 
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I. INTRODUCTION 



•i— > 

CO 

S ■ Over the years, polymers physics has greatly benefited from studies of lattice models. One persistent theme has been 
| the use of lattice models to uncover universal properties of chain molecules. An example is provided by the scaling 
exponents which characterize the statistical properties of polymer conformations, in the limit of very long chains [jp. 
For polymer chains confined to live in two dimensions, exact values of exponents were calculated by Nienhuis Q using 
the self-avoiding walk on the honeycomb lattice. The predicted value of the swelling exponent, which relates the linear 
size of the polymer to the number of monomers, was directly measured in recent fluorescence microscopy studies of 
DNA absorbed on a lipid bilayer . 
^ ■ Here we turn to the problem of polymer melting, which deals with a possible phase transition induced by the 
competition between chain entropy and bending rigidity. Bending rigidity determines the persistence length of the 
polymer. This is the distance over which the relative orientations of two chain segments are decorrelated due to 
thermal fluctuations. The long chain limit mentioned in the previous paragraph is obtained when the polymer length 
. is much greater than its persistence length. 
' It is important to point out that the effect of finite bending rigidity depends crucially on the steric constraints 
imposed on the polymer by its interactions with the solvent. For example, in the presence of a good solvent the 
polymer is in a "dilute" phase. Typical chain conformations are swollen with empty space between the monomers 
filled by solvent molecules. On the lattice, the dilute phase is characterized by a vanishing fraction of sites occupied 
H by monomers. In this phase, the bending rigidity simply increases the persistence length of the polymer, and it does 
not lead to a phase transition. This can be verified analytically in two dimensions, within the framework of Nienhuis' 
self-avoiding walk model § | . 

The picture changes considerably when the polymer is in a "compact" phase, with the monomers occupying all 
the available space. Such a situation is relevant, for instance, when modelling the conformations of globular proteins 
H . Compactness in this case follows from the interaction between hydrophobic amino-acids and the solvent (water) , 
. £h ! which leads to the expulsion of the solvent from the bulk of the protein. The simplest way to model this effect is to 
enforce compactness as a global, steric constraint on the polymer configurations pj. Within this compact phase, one 
expects a phase transition from a disordered melt to an ordered crystal as the stiffness of the polymer is increased. 

To study this melting transition, in 1956 Flory introduced a lattice model 0. Flory 's model, in its simplest 
formulation, consists of a single chain, described by a self-avoiding walk on the square lattice, endowed with a bending 
rigidity. To describe the melted phase the chain is taken to be maximally compact, filling all the sites of the square 
lattice; see Fig. jfl, The resistance to bending is modelled by an energy penalty for making 90° turns. 

In the Flory model, at infinite temperature the entropy dominates and the polymer will exhibit a finite density of 
bends, as in Fig. [l]a. As the temperature is lowered to zero all the bends are expelled from the bulk and their density 
goes to zero, as in Fig. The nature of the transition from the high temperature melt to the low temperature 
crystal has been debated over the years B. Here we show that the melting transition is continuous and calculate 
exact values of scaling exponents at the transition. 




FIG. 1: Compact polymer configurations on an 11 x 15 square lattice: a) Typical configuration in the melt phase, and b) 
zero-temperature crystalline state, in which the number of bends is minimum. 



In his original paper, Flory proposed a mean-field treatment which predicts a first order transition. According 
to J?]], the density of bends goes to zero at the transition and the chain entropy vanishes. This prediction of a first 
order transition with a vanishing entropy was challenged by Nagle Jo) . Namely, he showed that the exactly solvable 
six- vertex model maps to a related polymer model which differs from Flory's by the presence of polymer loops of all 
sizes. Applying Flory's mean-field approximation to this model leads once again to the prediction of a first order 
melting transition. However, as Nagle pointed out, this is at odds with the exact solution of the six-vertex model Jicf] 
which predicts a continuous, infinite order transition. This observation makes it questionable that the Flory approach 
is valid in the original model as well. In fact, a few years later Gujrati and Goldstein JHJ proved that the polymer 
entropy in Flory's model stays finite all the way down to zero temperature when it finally vanishes. However, the 
order of the transition still remained an unresolved question. 

Monte Carlo simulations of Baumgartner and Yoon |L2| , where they allowed for many chains and a finite density 
of empty sites, showed a first order melting transition. Soon thereafter Saleur |]l3| , using a transfer matrix approach, 
presented numerical evidence of a continuous transition, similar to the one found in the six-vertex model. More 
recently, Bascle, Garel and Orland (lj] proposed an improved mean-field treatment of the Flory model, which does 
not suffer from the problem of a vanishing entropy at the transition. It also predicts a first order transition. 

Here we show that polymer melting is continuous, as originally argued by Saleur fDjfl , by making use of a particular 
model, the semiflexible loop (SFL) model, and its height representation. Furthermore we calculate the central charge 
and exact scaling exponents at the transition. These results are checked against detailed numerical transfer matrix 
computations. 

The SFL loop model can be thought of as a "loop generalization" of the so-called F-model ||, in which suitably 
defined loops carry additional Boltzmann weights. The F-model is a special case of the six- vertex model |Tc|| , in which 
all vertices carry equal weights. This connection will serve as the motivation for introducing a more general model, 
the generalized six-vertex model, in which the general (zero-field) six-vertex model is endowed with extra loop weights. 
We shall finally introduce a similarly generalized version of the eight- vertex model |lC| ] . Its interest from a polymer 
point of view is that it allows for a unified description of semiflexible lattice polymers in a variety of phases: compact, 
dense and dilute. Furthermore it allows us to discuss the effect of vacancies on the polymer melting transition. 

The paper is organized as follows. In the next section we introduce the SFL model, which, in the limit of zero loop 
weight, gives the Flory model of polymer melting, and we discuss its phase diagram. In Sec. Ill we discuss the height 
representation of the loop model and how it leads to a conformal field theory in the scaling limit. We make use of 
the field theory in Sec. IV to calculate the central charge and scaling exponents, which we check against numerical 
transfer matrix computations in Sec. |v| In Sec. VI we propose a phase diagram for the generalized six-vertex and 



eight-vertex models. We end with a discussion of the scaling of semiflexible compact polymers, and we argue that 
the generalized eight-vertex model furnishes a rather complete description of non-compact semiflexible polymers. An 
appendix is reserved for a detailed discussion of the construction of the transfer matrices. 



II. SEMIFLEXIBLE LOOP MODEL 



Here we define the SFL model, and give a rough sketch of its phase diagram based on the limits of weak and strong 
bending rigidity. The fact that the SFL model reduces to the F-model in the limit of unit loop fugacity |J , plays an 
important role in guiding our intuition about the loop model. It also provides an exactly solvable line in the phase 
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diagram, against which the field theoretical and numerical results can be checked. 

A. Definition of the model 

The semiflexible fully packed loop model on the square lattice (the "semiflexible loop model" , or SFL for short) is 
defined by filling the square lattice with loops drawn along the lattice edges. Allowed loop configurations satisfy two 
constraints: 

• Self avoidance — loops are not allowed to cross, and 

• Full packing — every site is visited by exactly one loop. 

On the square lattice with periodic boundary conditions, edges that are not covered by loops also form loops, as there 
are two unoccupied edges associated with every site of the lattice. These we refer to as "ghost loops" . 

Given the configurations of the semiflexible loop model, the Boltzmann weights are defined in the following way. 
Every real loop is given weight rib, an d every ghost loop has weight n g . (In all the figures the real and ghost loops 
are shown as black and gray respectively, whence the subscripts b and g.) The parameters rib and n g act as fugacities 
of the two loop flavors, and as such they control the average number of loops of each flavor J3^| . They can be varied 
independently as the number of ghost loops is not fixed by the number of real loops Jl5| . Furthermore, a weight u>x 
is assigned to each vertex of the lattice at which the real and ghost loops cross. For wx > 1 this has the effect of 
disfavoring vertices at which the loop makes a 90° bend, or, in other words, the loops are semiflexible. The partition 
function of the semiflexible loop model is 

z = Y,< b < ew x, (1) 

Q 

where the sum runs over all allowed loop configurations Q. N\, and JV g are the number of real and ghost loops, 
respectively, while V is the number of crossing vertices; these are the two rightmost vertices in Fig. [|. In the limit 
rib —* 0, with rig = 1, we recover the Flory model: Z/n\, counts compact polymer loops each weighed by w^. 

The semiflexible loop model can be thought of as the generalization of the FPL 2 model introduced in Ref. [fL5| . 
The FPL 2 model is given by the partition function, Eq. ([!]), with wx = 1. It has a critical phase for |rib|, |n g | < 2, 
characterized by a power law distribution of loop sizes. For other values of the loop weights the model is non-critical 
with a distribution of loop sizes cut off at a finite value (fixed by the correlation length). Below we will show that 
the vertex weight wx, for each point in the critical phase of the FPL 2 model, produces a line of fixed points which 
terminates in a Kosterlitz-Thouless transition. 



B. Qualitative phase diagram 

Rough, qualitative features of the phase diagram of the semiflexible loop model can be deduced from the limits of 
zero and infinite bending rigidity. The motivation for developing a precise theory of the phase diagram, as mentioned 
in the introduction, stems from the interest in the rib ~ ► 0, n g = 1 case, which is the Flory model of polymer melting. 
We are also motivated by the relation of the SFL model to the integrable six-vertex model, and its generalizations. 

1. Flory model 

In the Flory limit of the SFL model, the wx = 1 point is the compact polymer problem, which we have studied 
previously Jig] , Here one is concerned with enumerating all self-avoiding walks that visit every site of the lattice. 
We have shown that compact polymers on the square lattice are a critical geometry characterized by non-mean-ficld 
scaling exponents which can be calculated exactly from a field theory. 

As u>x is increased away from one, we are dealing with a compact polymer with a bending rigidity. In the limit 
wx — > oo we arrive at a frozen phase in which the density of vertices at which the polymer bends goes to zero. This 
is the polymer crystal. At an intermediate weight wx = w x (1 < u> x < oo) there will be a melting transition. One of 
the important unresolved problems is the nature of this transition. Here we construct an effective field theory of the 
Flory model and show that the melting transition is continuous. 

Another interesting issue is the region of < wx < 1 . As wx — * straight-going vertices are completely suppressed, 
and with appropriate boundary conditions the only allowed configurations are those of a checkerboard pattern of small 



FIG. 2: Typical configurations in the SFL model with rib = n s = 2 and bending rigidity parameter wx = 1/4 (left panel), 
Wx = 1 (middle) and tux = 4 (right). We shall show that the left and middle panels correspond to critical melt st ates, w hile 
the right panel is a non-critical crystalline state. In the latter, domains of non-zero staggered polarization (see Sec. II B 2| ) are 
clearly visible. 



FIG. 3: Correspondence between the vertices of the six-vertex model and the FPL 2 model (here shown for an even vertex; at 
odd vertices the arrows are reversed). 

loops, each loop having its minimal length of four. If the Flory limit (rib — > 0) is taken before the wx — > limit, there 
has to be a number of straight-going vertices at the boundary, the dominant configurations being those of a single 
wiggly line. In any case, the wx ^ limit is again a crystalline phase of zero entropy. We shall however argue below 
that the corresponding crystallization transition is located at wx — and is thus rather uninteresting. 

A qualitative idea of the physics underlying the phase diagram of the SFL model can be obtained by looking at 
some typical configurations for various values of wx', see Fig. ||. The images were obtained by performing Monte Carlo 
simulations on a square lattice of size 100 x 100 with toroidal boundary conditions. For technical reasons [Q we take 
rib = n s — 2 and no loops of non-contractible topology are allowed. (Further details on the algorithm used for these 
simulations can be found in [fL6|.) 



2. Six-vertex model 



Before turning our sights to the semiflexible loop model it is instructive to review exact results for the (zero-field) 
six- vertex (6V) model. The 6V model corresponds to the = n g — 1 line in the phase diagram of the SFL model. 
The mapping between the two is simple: at even (odd) vertices the edges covered by the real loops are identified with 
arrows pointing out (in), while the edges covered by the ghost loops correspond to arrows pointing in (out); see Fig. ||. 
The appropriate six- vertex weights are: a = b = 1 and c = wx [ f33| . 

In the 6V model there is an order-disorder transition as a function of the vertex weight wx- In the ordered state, 
which is obtained for wx — > oo, all the vertices are of the c- variety (cf. Fig. ||). The order parameter is the staggered 
polarization, which in the loop language can be expressed as the difference between the number of horizontal and the 
number of vertical loop-covered edges, per site Ol. The exact solution of the six-vertex model predicts a continuous 
(infinite-order) transition occurring at wx = 2 ||l0|. The disordered phase for < wx < 2 is critical with an 
infinite correlation length and power-law correlations. Below we will show that there is an analogous transition in the 
semiflexible loop model, as wx is varied, for all values of |rib|, |n g | < 2, including the Flory case (rib - * 0,n g = 1). 
In the Flory model this was observed previously by Saleur in numerical transfer matrix computations p3[. For the 
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critical phase of the model we shall construct an effective field theory using the interface representation of the loop 
model. This leads to exact (but non-rigorous) results for the central charge and scaling dimensions, which we confirm 
via numerical transfer matrix calculations. 

III. FIELD THEORY CONSTRUCTION 

To construct a field theory for the critical phase of the semiflexible loop model (SFL) we make use of the height 
representation of the fully packed loop model on the square lattice (FPL 2 ). This was already described in detail in 
our previous work fl5|| , and here it is briefly reviewed for completeness. The main effect of the vertex weight wx on 
the field theory is to renormalize one of its coupling (elastic) constants. This does not change the central charge, but 
leads to continuously varying scaling dimensions for a specific subset of operators, which we identify. These results 
are confirmed by our numerical transfer matrix computations. 

A. Height map 

The height mapping is defined on the space of oriented loop configurations {G 1 }- We associate 2 Ne+Nh oriented 
loop configuration Q' with each loop configuration Q of the SFL model by independently orienting every real and 
every ghost loop clockwise or counterclockwise. 

The Boltzmann weight of an oriented loop is exp(i0), where the phase 4> = ±7reb for clockwise (counterclockwise) 
oriented real loops, and <j) — ±7re g for the two orientations of the ghost loops. To recover rib and n g for the loop 
weights, after summing over the two possible orientations we must set 

rib = 2cos(7reb) 

n g = 2cos(7re 9 ) . (2) 

This particular partition of the loop weights between the two orientations has the advantage of allowing the loop 
weights to be distributed among all the vertices that the loop visits, thus rendering the weights local. This is achieved 
by assigning the phase 7reb/4 (7re g /4) to every vertex at which the oriented real (ghost) loop makes a right turn, and 
the opposite phase for left turns. The fact that for every closed loop on the square lattice the difference between the 
number of left and right turns is ±4, is what makes these vertex weights work. The total vertex weight is then given 
by the product of phase factors when the loops bend, while a weight % is assigned to vertices at which the loops do 
not bend. 

Turning back to the height map, we define microscopic heights h(x) on the lattice {x} dual to the square lattice 
on which the loops are defined. Once the height at the origin is fixed, the heights on all the other vertices of the 
dual lattice are uniquely specified by the oriented loop configuration. Namely, the height difference between nearest 
neighbor vertices of the dual lattice is A, B, C or D, depending on the state of the edge that separates them. The 
four height-difference vectors, also referred to as "colors" , are associated with the four possible states of any given 
edge, which can be either covered by a real or a ghost loop, with one of two possible loop orientations. Real loops 
are formed by alternating cycles of A and B colored edges, while the C and D colored edges are ones visited by the 
ghost loops. Note that the difference between an ABAB • • • and a BABA • • • cycle encodes the orientation of the 
corresponding (real) loop. 

The fully packing constraint and the requirement that the height be unique (i.e., the sum of height differences along 
any closed lattice path must be zero), imposes a single algebraic constraint on the four colors: A + B + C + D = 0. ft 
follows that only three of the four vectors are linearly independent. A convenient choice that respects the symmetries 
between the four colors is to let the corresponding vectors point from the center to the vertices of a regular tetrahedron: 

A = (-1, +1, +1) , B = (+1, +1, -1) , 

C = (-1,-1,-1) , D = (+1,-1, +1) . (3) 

The effective field theory for the SFL model describes the fluctuations of the coarse-grained heights which retains 
only the long- wavelength (much larger than the lattice spacing) Fourier modes of the microscopic heights. 

B. Effective field theory: i»x = 1 

For wx = 1 we have the familiar case of the fully packed loop model on the square lattice. Its effective field theory 
was discusses in a previous publication Jig] and here it is reviewed for completeness. 
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The partition function of the loop model in the height representation can be written as a path integral over the 
coarse grained heights with the (dimcnsionless) action: 

S = Se + S b + S l . (4) 

This action only takes into account the long- wavelength fluctuations of the microscopic height. The three terms in 
the action are of different origin. 
The elastic term, 

S E = \ J d 2 x {Knftdh 1 ) 2 + (dh 3 ) 2 ] + 2K 13 (dh 1 ■ dh 3 ) + K 22 {dh 2 ) 2 } , (5) 

accounts for the height fluctuations due to the entropy of fully packing the square lattice with oriented loops. Equiva- 
lently, this is the entropy of edge coloring the square lattice with four different colors. The elastic term favors oriented 
loop configurations that minimize the variance of the microscopic height; these are the macroscopically flat states. In 
terms of the color degrees of freedom the fiat states have the property that the four edges of each elementary plaquette 
are colored by two colors only. 

The particular form of the matrix of elastic constants, K, is fixed by the lattice symmetries and symmetries 
associated with permuting the colors A, B, C, and D. The elastic constants Kg are functions of the loop fugacity. 
For the wx = 1 case they were calculated in Ref. [[l5| using the loop ansatz |l7j], which allows one to identify the 
marginal screening charges [pTsfl . For the FPL 2 model there are four screening charges: 



(!) = (-^0, +7T) 



e 

= (-7T,0,-7r) 

e (3) = (-7T, +7T,0) 

e< 4 ) = (-7r,-7r,0) . (6) 

These electric charges are associated with the most relevant vertex operators appearing in the Fourier expansion of 
the operator conjugate to the loop weight (see Eq. (|l2]), below). 

Demanding that all four charges have scaling dimension equal to two gives (using Eq. |l6| ) 

Ku = g (2 - e g - e b ) 

Kis = l(e b ~e g ) (7) 
7r (l-e 6 )(l-e g ) 

A 22 — 

2 2 — e b — e g 

for the elastic constants of the FPL 2 model; eb and e g satisfy Eq. (||) and take their values on the interval [0, 1/2]. 
Below we will argue that the effect of wx 7^ 1 is to change the value of the elastic constant K22 while leaving the 
other two unchanged. 

The boundary term in the action, 



l - J d 2 x (e • h(x))p(x) , (8) 



5 B = — 
Air 

enforces the correct weight of topologically non-trivial loops. If the oriented loop model is defined with periodic 
boundary conditions along one direction (i.e., on a cylinder) these would be the loops that completely wind around 
the cylinder p4| . On a cylinder the scalar curvature, p, is non-zero only at the two boundaries at infinity. Sb has the 
effect of placing background electric charges ±eo at the two boundaries, where the identification 

7T 

e = --(e g + e b ,0,e g -e b ) (9) 

comes about by demanding that the oriented winding loops be assigned correct phase factors, exp(±i7reb) or 
exp(±i7re g ) |§. 

The third term, called the Liouville term, 

S L = /d 2 x w[h(x)] , (10) 
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owes its existence to the c omple x weights associated with oriented loops in the bulk. The local redistribution of the 



loop weights made in Sec. Ill A leads to complex vertex weights, which in turn depend only on the colors of the four 



edges around the vertex. If we write the vertex weight as exp(— w) then 

iu(B,C,A,D) = 0, 
iu(B,D,A,C) = 0, 



. 7T 



w(A,B,C,D) = T i-(e g + e b ) 
w(B,A,C,D) = =pi^(e g -e b ) 



10(A,B,D,C) = Ti^(e b -e g ), 

«;(B,A,D,C) = TiT(-e b -e g ); (11) 



the top (bottom) sign is for even (odd) vertices, and the colors are listed in order, starting from the left-most edge 
and proceeding clockwise around the vertex. The weight operator w is invariant under cyclic permutations of the 
colors and it is a periodic function of the heights around a vertex. In the scaling limit the vertex weights give rise to 
the operator w[h(x)] in Eq. ( |To| ) which can be written as a Fourier series 

Mh(x)] = ^2 * e exp(ie • h(x)) . (12) 

The electric charges e appearing in the Fourier expansion are dictated by the lattice of periodicities 1Z W of the operator 
w[h]; IZ^ is the reciprocal lattice. 1Z W is determined by inspection of the values the loop weight operator takes on 
the flat states: vectors in 1Z W connect flat states on which the loop weight operator takes identical values. The most 
relevant charges in Ti* w are the four given in Eq. @ . We identi fy t hem with the screening charges [Q of the Coulomb 
gas. This is the content of the loop ansatz introduced in Rcf. |17j| . 



C. Effective field theory: wx 7^ 1 

For the SFL model, when wx ^ 1, the Liouville term in Eq. (0) is modified, while the elastic and the boundary 
terms are unchanged. The number of marginal screening charges appearing in Eq. (|l2| ) is reduced from four to two, 
and the loop ansatz fixes the values of K\\ and only. They do not depend on the value of wx and are given by 
the wx = 1 formulae, Eq. (Q). K 22, on the other hand, is a non-universal function of wx- Below we present arguments 
for this scenario, which is supported by exact results available in the 6V case (i.e., for n b — n g — 1) 5 an d by our 
numerical transfer matrix calculations described in Sec. |v| below. 

The new vertex weight wx changes the value of w in Eq. ( pi] ) from to — Inwx for the vertex states (B, C, A, D), 
(B,D, A, C), and six other related to these two by cyclic permutations of the colors. The weights of the other 16 
vertex states are unchanged. We consider the consequences of this change on the effective field theory. 

In the height representation of the SFL model, the change in vertex weight corresponds to adding 

Sx = y"d 2 x X[h(x)] , (13) 

to the action. The X operator takes the value In wx on the flat states made up of (B, C, A, D) or (B, D, A, C) type 
vertices, and vanishes on all the others. By inspection of the graph of flat states we find that the lattice of periodicities 
for the operator X, Tlx, is the span of (1, 0, —1), (1, 0, 1) and (0, 1, 0); these are the height difference vectors between 
the flat states in the support of X. This observation implies that X[h] can be expanded in a Fourier series over 
electric charges that live in the dual lattice TZ X which is the span of (it, 0, ir), (tt, 0, — it) and (0, 2ir, 0). 

If we consider the effect of Sx as a perturbation on the action of the FPL 2 model the electric charges (0, ±2ir, 0) 
play a special role. Namely, the operator product expansion of exp (i(0, 2n, 0) • h) and exp (— i(0, 2tt, 0) • h) contains 
the (dh 2 ) 2 operator, and therefore leads to the renormalization of K22 pl[ . This follows from the fact that the 
background charge, Eq. (^|), has a vanishing second component. On the other hand, for charges e with non-zero first 
or third component, the effect of the background charge is that the operator product expansion of exp (ie • h) and 
exp (— ie ■ h) does not contain dh 1 ■ dh? operators and therefore does not lead to the renormalization of the clastic 
constants K*j. 
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The Coulomb gas representation of the height model provides a clear physical picture of the effect of Sx on the 
critical action of the FPL 2 model. For wx — 1 the dimension of the (0, ±2w, 0) charges follows from Eq. (|l6|), 



471^22 \ 1 — e b 1 

It is greater than 2 in the whole critical region of the FPL 2 model. These charges are therefore irrelevant in the 
renormalization group sense. In the Coulomb gas picture the (0, ±2tt, 0) charges appear as bound pairs of neutral 
dipoles. Increasing wx will have the effect of increasing the bare fugacity of these dipoles which will in turn increase 
the value of the coupling K22 appearing in the effective field theory. Formally, this can be seen in perturbation theory 
making use of the operator product expansion pl| . Physically, the renormalization of if 22 can be understood as the 
screening effect of dipoles. The dipoles lower the Coulomb energy between two electric test-charges having a non-zero 
second component, corresponding to an increase in the value of K22 which plays the role of a dielectric constant. 

At a critical value w x there will be a Kosterlitz-Thouless type transition of the SFL model into a flat state with a 
vanishing density of vertices at which the polymer bends. At the transition the (0, ±27r, 0) charges are marginal, i.e., 
their scaling dimension is equal to 2. Using Eq. (Q) this observation gives rise to the prediction for the critical value 

of K 2 2- 

K22(W X ) = \ . (15) 

For values of wx smaller than io|, if 22 will be a non-universal function of wx- In the nb = n g = 1 case, the formula 
if 22 = arcsin(u>x/2) follows from the exact solution of the 6V model |L0]]. The critical value of the vertex weight is 
w x (6Y) — 2 and if 22 (2) = tt/2 is in agreement with Eq. ([l5|). For other values of nt, and n g our numerical transfer 
matrix calculations are in good agreement with Eq. (|l5|). 

The introduction of the vertex weight wx also has an effect on the screening charges, Eq. (^|), that appear in the 
Liouville part of the action. First consider the rib = n g case of the SFL model. Due to the presence of the wx term 
cyclic permutations of the four colors around a vertex are no longer a symmetry of the vertex weight. Therefore, 
unlike the wx — 1 case [jlBI , there are now two independent elastic constants, if 22 an d ifu appearing in Sb- if 13 = 
follows from the remaining Z2 symmetry of the vertex weights which are invariant under two cyclic permutations, 
such as (A, B,C,D) — * (C,D,A,B). The deduced structure of the elasticity matrix implies that the four electric 
charges in Eq. (^|) are no longer degenerate in dimension for arbitrary wx- Since the dimensions of e^ 1 ) and e*- 2 -* 
are independent of ^22 they are identified as the two screening charges tied to the non-renormalizability of the loop 
weights [ p5| . As in the FPL 2 model we then assume that these two charges remain marginal when 7^ n g . Using 
the dimension formula, Eq. (|l6|), this then fixes the values of the two elastic constants, K\\ and if 13, to the values 
quoted in Eq. (|7|). 

Finally, it is interesting to look at some extreme limits of if 22 in view of the effective field theory. Consider first 
the limit if 22 — > 00 hi which height fluctuations in the second height component are completely suppressed. (As we 
are outside the critical phase, we are here referring to the bare value of the coupling.) Clearly, height fluctuations 
must always be present in the microscopic four-coloring model, but it is nevertheless instructive to look for the states 
that minimize the fluctuations of h 2 . From the choice of the color vectors, Eq. (|^), it is not difficult to see that on 
the four sites of {x} surrounding a given vertex, h 2 fluctuates by two units for the first four vertices of Fig. || and by 
one unit for the last two vertices. All vertices must therefore be of the c-type, corresponding to the limit wx ~^ 00. 
Thus, K22 — * 00 as Wx — * 00. 

Conversely, as if 22 — > the fluctuations in h 2 become unbounded and the effective field theory loses its consistency 
(since it was based on the assumption that the interfacial entropy is due to bounded fluctuations around the macro- 
scopically flat states). However, the argument given above indicates that a small value of K22 should correspond to a 
small number of straight-going vertices in the loop model. Thus, we would conjecture that if 22 — * as wx — > 0. This 
expectation is confirmed by the exact result for the 6V case |10[ and also by extrapolation of our numerical results 
for K22{wx) in the Flory case. 

Apart from these limiting values, we would of course expect if 22 to be a monotonically increasing function of wx 
throughout the critical phase. 

In the next section we compute the central charge and the scaling dimensions of various operators in the semiflexible 
loop model from its effective field theory. We identify quantities that depend on the non-universal clastic constant 
if 22; these are then predicted to vary continuously with wx- 



IV. OPERATORS AND SCALING DIMENSIONS 



The effective field theory of the semiflexible loop model describes a Coulomb gas of electric and magnetic charges 
in the presence of background and screening charges. The magnetic charges m are vectors in 1Z which is the lattice 
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of periodicities of the graph of flat states, while the electric charges e take their values in the reciprocal lattice TV 
p5|. With the normalization adopted for the vectors A through D, Eq. 72. is a face-centered cubic lattice whose 
conventional cubic cell has sides of length 4, while is 1Z* is a body-centered cubic lattice whose conventional cubic cell 
has sides of length ir. 

The scaling dimension of an operator which has total electromagnetic charge (e, m) is the sum of its electric and 
magnetic dimensions, and it is a function of the elastic constants and the background charge |L8| : 

x(e,m) = — [(eKT 1 ) . ( e - 2e ) + (mK) • ml . (16) 

4-7T L J 

K is the 3x3 matrix of elastic constants and K" 1 is its inverse. 

From Eq. @ and the form of K (Eq. (§)) and eo (Eq. (0)) it immediately follows that operators whose electric 
and magnetic charges both have a vanishing second component will have a /^-independent scaling dimension. The 
scaling dimension in this case is independent of w% and equal to its known value at wx = 1 JT^ . Operators with e 
and m charges whose second components are not both zero will, on the other hand, have a scaling dimension that 
varies continuously with wx- These predictions are confirmed by our numerical results. 



A. Central charge 

The central charge of the SFL model follows from its critical action. The three height components (bosonic free 
fields) each contribute one unit to the central charge while the contribution from the background charge is 12x(eo, 0). 
Using Eq. ( fL6| ) for x(eo,0), Eq. (||) for the background charge eo, and the calculated values of the elastic constants 
Kn and K 13 , Eq. (0), we find 

c = 3-6f-^- + -^V (17) 
y 1 - e b 1 - e g J 

independent of the unknown value of K 22 ■ 

For the 6V model, which corresponds to the eb = e g = 1/3 line in the SFL model, the above formula gives c = 1 
for the central charge along the critical line. This result also follows directly from the exact solution of the 6V model. 

For the Flory model of polymer melting, which is the eb = 1/2, e g = 1/3 case, thepredicted central charge is 
c = — 1. This value is confirmed by our numerical transfer matrix calculations (see Sec. [VJ). 



B. Thermal operator 

The SFL model can be thought of as the zero-temperature limit of a more general model where we allow for thermal 
excitations that violate the fully packing constraint. Violations of the constraint lead to vertices with the four adjacent 
edges colored (C,D,C,D). In the height representation such a vertex is identified with a topological defect (screw 
dislocation) whose charge, i.e., the sum of height differences around the vertex, is 

m r = 2(C + D) = (0,-4,0). (18) 

Other vertices which have no A or B colored edges are possible, but they have a larger magnetic charge and are hence 
less relevant. 

In the Coulomb gas picture a topological defect corresponds to a magnetic charge. Therefore, the thermal dimension 
can be calculated using Eq. @, and we find 

x T = z(0, m T ) = -K 22 ■ (19) 

7T 



We make use of this equation below as it allows us to determine the unknown elastic constant K 22 from a measurement 
of the thermal scaling dimension. Once this elastic constant is known, scaling dimensions of all electromagnetic 
operators can be calculated from Eq. (16). 



10 



C. String operators 



A particularly important set of operators in any loop model are the string operators. Their two-point function 
is defined as the probability of having the small neighborhoods around two fixed points on the lattice, which are 
separated by a large distance, connected by Sb real loop segments and s g ghost loop segments. For simplicity, we 
shall require Sb and s g to be either both even or both odd; Sb + s g odd requires L to be odd which produces a twist 
in the height, as discussed in Ref. fl5|| . In the height representation these string configurations are mapped to two 
topological defects, one serving as the source and the other as the sink of oriented loop segments. When the oriented 
loop segments wind around the defect points they are assigned spurious phase factors by the vertex weights; these 
pha se factors can however be compensated by introducing appropriate electric charges at the positions of the defects 

In the case Sb = 2&b and s g = 2fc g , i.e., when the number of real and ghost strings are both even, the electric and 
magnetic charge of the corresponding string operator are |i~5f| 



e2fc b ,2fe g 



2 

-2(fe b 



(e b ,0, -e b )(l - 4 b .o) - -(e g ,0,e g )(l - 5 k ) 



(20) 



Since the charges have vanishing second component their dimension is independent of K22 and constant along the 
whole critical line wx < w x . The value of the string dimension follows from Eq. (|l~6|), 



X2k b .2k„ 



(1 - e h )ki + (1 - e g )kl - 



1 - e b 



(l-4 b ,o)- 



1-e, 



■(l-4 g ,o) 



(21) 



and is identical to that of the FPL 2 model ]l~5[ . Our numerical simulations confirm that even string dimensions are 
constant along the critical line. 

In the odd string case, when s b = 2fcb — 1 and s g = 2fc g — 1, the electric and magnetic charge are |fl5f 



7f 

e 2 fc b -i,2fc g -i = -^{e g + e h ,0,e g - e h ) 
m 2 fc b _i,2fc g -i = -2(&b + k s - 1, 1, feg - kb) ■ 
Notably the magnetic charge has a non- vanishing second component. Using Eq. (n&j we calculate 



(22) 



z 2fcb -i,2fc g -i = ^ + I [(1 - e h )(2k b - l) 2 + (1 - e g )(2fc g - l) 2 ] - \ 



(23) 



for the odd string dimension. It depends on the value of ^22 and will therefore vary continuously with wx- At the 
melting transition the exponents are exactly known from Eq. ( |l5[ ) . This is confirmed by our numerical transfer matrix 
results, which we describe next. 



V. TRANSFER MATRIX RESULTS 



To check the correctness of our field theoretical predictions, we have numerically diagonalized the transfer matrix 
of the semi-flexible loop model (and of its various generalizations, to be discussed below) defined on semi-infinite 
cylinders of even widths L ranging from 4 to 14. 

The existence of a transfer matrix may not be a priori obvious, since the Boltzmann weights depend on the number of 
loops, which is a non-local quantity. We have however already shown in an earlier publication p"5| how this is resolved 
by working in a basis of states that contains non-local information about how loop segments are interconnected at 
a given stage of the computations. In that paper, it was also shown that the full transfer matrix contains various 
sectors, the leading eigenvalues of which provide finite-size estimations of the free energy and of the various critical 
exponents, using the standard CFT relations [^4|, |25| 

f (L) = / (oo)-^+-. (24) 

ML) -f (L) = (25) 
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FIG. 4: Central charge c as a function of the bending rigidity wx for the 6V model (left) and the semiflexible loop model 
(right). We show three-point fits for different system sizes, as indicated on the legends. 



Here, the label k refers either to a higher eigenvalue in the sector to which /o belongs, or to the leading eigenvalue in 
another sector characterized by some topological defect of charge k. 

To access critical exponents, we shall mainly be concerned with topological defects that consist in enforcing that 
a certain number of strings of either flavor propagate along the length direction of the cylinder. These give rise 
to a two-parameter family of critical exponents x Sb ^ s corresponding to Sb real strings and s g ghost strings. The 
corresponding topological charges, Eq. (^0|) and Eq. (^2|), take the form of three-dimensional electromagnetic vector 
charges. In the transfer matrix calculations, each of these topological sectors is associated with a different state space. 
The difficulty of precisely characterizing these spaces limited our previous approach [|l5| to at most two strings. In 
Appendix |a| we present an algorithm that explicitly constructs the required state spaces for any (sb, s g ), based on an 
iterative procedure and hashing techniques. 

By inspection of the eigenstates produced by our previous algorithm [jl5| , it turns out that many of the basis states 
carry zero weight. One would then expect that identical results can be obtained more efficiently by working in a 
basis in which such states have been eliminated from the outset. We defer the technical details of how this can be 
done to Appendix |X[ It is also shown how the block-diagonalization scheme can be carried even further, by exploiting 
various conservation laws that are most easily understood from the analogy between the SFL model and the six- vertex 
model. One important consequence is that the constrained free energy fr(L) that is linked to the thermal scaling 
dimension can now be obtained as a leading eigenvalue, rather than as the second eigenvalue in the stringless sector. 
This considerably improves the efficiency of the computations. 

Finally, the matrix elements need some modification in order to take into account the bending rigidity parameter 
wx- This is readily done, without any modification of the basis states, since wx is a purely local quantity. 

Before turning to our numerical results, we should mention that we have submitted our transfer matrices to several 
tests, in order to verify their correctness: 

• For wx = 1, all numerically determined string dimensions x Sb _ Sg with Sb + s g = 2 or 4 agree to at least three 
significant digits with their exact values in the cases (rib, n g) = (lj 1) 10 an d («b, %) = (0, 1) |lq| . 

• All eigenvalues found for the FPL 2 model agree with those obtained from our previous algorithm |l5| . 

• For the six- vertex model JlO[ , we have compared the extrapolated bulk free energy with Baxter's exact expression. 

• Again for the six- vertex model, we find excellent agreement with the exact formulae x\^\ — -^p- and xt = 4J ^ 22 , 
where K22 = arcsin(ux/2) is the elastic constant. 

• We have also found agreement with the first few terms in diagrammatic expansions around various limits of 
infinite fugacities. 



A. Central charge 

A crucial prediction of our field theory is that, for given values of the loop fugacities (rib, rig), the central charge of 
the SFL model should be independent of wx, as long as the latter is constrained to the critical regime, < wx < w x . 

In Fig. H we show the effective central charge as a function of in the cases (n-^,n g ) — (1, 1) and (0, 1). The 
result c — 1 for < wx < 2 is well established for the 6V model, but the plot for this case is still useful as it gives 
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us some guidance as to what finite-size effects to expect. In particular, note that these become more pronounced 
when u>x is small, and that the termination of critical behavior at wx — 2 is clearly signalled by the finite-size data's 
levelling off. Another effect is that while for wx < 2 the distance between successive estimates decreases with system 
size, for wx > 2 we observe this distance to increase. 

Although finite-size effects play a more important role in the (nb,n. g ) = (0,1) case, the general picture is quite 
similar. The figure leaves little doubt that c = — 1 for < wx < w x . We also obtain a first rough estimate 
w x = 1.95 ± 0.15, not far away from the 6V-model value. 

Here, and elsewhere, we mainly show fits in which the convergence of the relations Eqs. ( gjj) and (p5|) has been 
accelerated through the inclusion of a non-universal 1/L 4 correction, as predicted by conformal invariance. 



B. String dimensions 

We next turn to the computation of the magnetic- type scaling dimensions x Sb . s describing the scaling of the 
operator that inserts Sb real strings and s g ghost strings. To study these, the width L of the strip must have the 
same parity as Sb + s g . For simplicity we shall limit ourselves to the case of even L. There are then two classes of 
exponents: Those in which Sb and s g are both even, and those in which they are both odd. The field theory predicts 
that the former should stay constant on the critical line, parameterized by wx, while the latter are expected to vary 



continuously as functions of wx] see Sec. [VC 

On Fig. I we show two examples of exponents with Sb,s g even, within the SFL model [(rib,%) = (0,1)] with 
varying wx- They correspond respectively to the insertion of two ghost strings (2-0,2) and of four real strings [x^q). 
From the figure it should be evident that £0,2 = \ an d 2:4,0 = f ar e constant throughout the critical phase. In the 
latter case the finite-size variations are quite pronounced, as might have been anticipated given the higher number of 
strings. Careful observation of the distance between subsequent finite-size points however strongly suggests that the 
variation will eventually die away. 

Examples of exponents with Sb,s g odd are given on Fig. g. In both cases, 2,1,1 and 2:3,1, the convergence to 
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TABLE I: Thermal exponent xt measured for varying values of wx- The corresponding values of the scaling dimensions xi,i, 
£3,1 and 2^3 are compared with their directly measured counterparts. 



monotonically increasing functions of wx is clearly brought out. Also note the agreement with the exact results for 
wx — 1, which read respectively x\ t \ = — yfj — — 0.0446 and x 3 ,i — — 0.455 Q |. 

C. Thermal scaling dimension 

As described in Appendix |A|, the thermal scaling dimension is linked to the gap between transfer matrix sectors in 
which there is an even (respectively an odd) number of flavor crossings in the basis states. Because of the relation 

K 22 = Jx T , (26) 

measuring this gap gives a direct means of accessing the elastic constant associated with the second height component 
in the field theory. 

According to the field theory, K22 is a non- universal function of wx, and once it is known the values of all the other 
critical exponents follow. This suggests the following numerical check of the field theoretic scenario: For several values 
of wx, we measure xt from the transfer matrix, and use it to determine K22- We then compute the predictions for 
the various other scaling dimensions (the %s b ,s s ) from the field theory, Eq. ([l6|), by use of the numerically determined 
value of K22, and compare them with values measured directly from the transfer matrices. 

The result of this verification is shown in Table || The values for xt are based on transfer matrices for strips up 
to size L = 14, here extrapolated to the limit L — > 00. The agreement between the CFT predictions and numerics 
is in general excellent. Note however, that the precision deteriorates whenever wx approaches zero or w x , its critical 
values. 

Based on these data we can refine our estimate for the location of the melting transition: 

w x = 1.92 ± 0.02. (27) 

In Fig. ^ we compare our numerical results for the curve xt{wx) in the SFL case with the exactly known result of 
the 6V-model, xt = ^ arcsin(u>x/2) |lC| ]. Although the functional forms are quite reminiscent, we have unfortunately 
not been able to conjecture a convincing exact expression in the SFL case. 



14 



2.0 F 




0.0 0.5 1.0 1.5 2.0 



FIG. 7: Thermal scaling dimension xt versus bending rigidity wx in the Flory model (symbols), as compared to the exact 
result of the six- vertex model (line). 



VI. PHASE DIAGRAM 



A. Generalized six-vertex model 



Given the one-to-one correspondence between the six vertex configurations in the FPL 2 model and the six arrow 
configurations in the six- vertex model (see Fig. |^), it is natural to define a generalized six- vertex model in which the 
standard arrow weights are supplemented by the non-local loop weights rib, n g of the FPL 2 model. 

Until now we have only considered the isotropic case of a — b (see Fig. Q. Let us briefly recall the effect of taking 
a b in the six- vertex model JlO[ . Define the parameters w and /i by 

a 2 + b 2 — w'i „ 

A = — ^ = -cos/i, 0</i<7r 28) 

2ab 

a exp(i/j) - exp(iw) 

T = / ■ ~T" ■ \ — T' — ^ < w < (i. (29) 

o exp(z/i + iw) — 1 

Then, taking a ^ b corresponds to twisting the usual square lattice into a rhombus, defined by the anisotropy angle 

n 

6 = -{! + -). 30 
2 /i 

All this means is that the central charge and the critical exponents, when measured in the usual way from a transfer 
matrix, get multiplied by a geometrical factor of sin(#). 

In Fig. H we plot the effective central charge of the SFL model with b = wx = 1 and varying a against the variable 
t = 6/tt, defined in terms of the above 6V expressions. By the word "effective" we mean that we do not correct for 
the lattice distortion, the effect of which can then be read off from the graph. If the effect of the anisotropy were the 
same as in the 6V model, the plot should just look like the function — sin(7rr), since the SFL model has (real) central 
charge c = — 1. Clearly, this is not the case, and so the non-locality of the loop weights has a non-trivial effect on the 
anisotropy factor. We leave this as an interesting open question. 



B. Generalized eight- vertex model 

It is also of interest to consider the loop generalization of the eight-vertex model. In terms of the loops there are 
two different ways of resolving the vertices that act as sources or sinks of the eight-vertex arrows, and so we are led 
to consider the ten-vertex model defined by Fig. ^. In addition to the local weights which are shown on the figure, 
we assign the usual non-local loop weights rib and n g . 

For simplicity, we shall disregard the effects of anisotropy, and thus only two types of local weights are of interest. 
The first is the weight wx of having the two loop flavors cross, same as in the SFL model. The second is a contact 
interaction w c , assigned to the vertices where two loop segments of the same flavor touch one another. One may 
consider letting it depend on the flavor index, but in order to stay close to the definition of the conventional eight- 
vertex model we shall here take the contact interaction to be flavor independent. 
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FIG. 9: Vertices defining the generalized eight-vertex model, along with their corresponding multiplicities and local weights. 



The motivation for the contact interaction is to be able to exclude the loops of a given flavor from any number of 
lattice vertices. As this violates the compactness constraint, we expect the conclusions of our earlier paper on the 
transition from the compact to the dense phase [p3| to apply: A non-zero value of w c should induce a flow towards 
a phase where the two loop flavors decouple, and the critical properties are just those of two non- interacting 0(n) 
models (with n — and n g respectively) in the low-temperature (dense) phase. 

A detailed numerical study of the behavior of the effective central charge in the parameter space (wx, w c ) has led 
us to suggest that the phase diagram of the generalized eight-vertex model is as shown on Fig. lp| . 

For w c = 0, the model reduces to the SFL model, and so below the melting point M (i.e., for wx < w x with, very 
roughly, w| ss 2) we have a line of critical points along which critical exponents that depend on the second height 
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FIG. 10: Proposed phase diagram of the generalized eight-vertex model. 
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component vary continuously, while the central charge 

6el 6e 2 

cK ng )=3--^--A- (31) 
1 - e b 1 — e g 

is constant. The end-point of the SFL line, with wx = oo, is a trivial attractive fixed point T, favoring configurations 
in which all loops go straight in the bulk (they are necessarily reflected at the boundaries enjoying free boundary 
conditions). The point T is believed to govern the theories to the right of M (i.e., with wx > w x ), including a portion 
of the phase diagram with non-zero but small w c (see Fig. |l0| ). 

Moving away from the critical line of the SFL model, towards positive values of the contact interaction, we observe 
numerically that the central charge drops abruptly by one unit, and stays constant as a function of w c up to some 
finite critical value u£ that depends on wx- This is the dense phase of the DPL 2 model |^3| with central charge 

C( „ b ,„ g)= ( 1 -_M_) + ( 1 - T M_). (32 ) 

Here, the two loop flavors decouple, and critical exponents are just the sum of the critical exponents for two non- 
interacting O(n) models (with n — rib and n g ) in the dense phase. We have verified numerically this prediction for the 
exponent x\ t \ for a number of different loop fugacities. We have also observed numerically that the critical exponents 
do not depend on wx throughout the dense phase. This confirms the expectations that in non-compact phases the 
only effect of the bending rigidity is to rcnormalize the persistence length of the polymer, as already discussed in the 
introduction. 

Finally, for w c large enough, the numerically evaluated central charges suggest that the models flow into an attractive 
fixed point P situated at (wx,w c ) = (0, oo). Here, only contact-type vertices are allowed, and since the different loop 
flavors can no longer coexist, the partition function at P becomes simply a sum, Z = Zb + Z s , where Zk involves only 
contact vertices of flavor k (with k = b,g). But clearly Zk is just the loop- model representation (l(| of a self-dual 
Potts model with qk = (rik) 2 states. It is intuitively clear (and explicitly brought out by the exact solution ]l0|]) that 
the free energy of the g-state Potts model is an increasing function of q. Therefore, the sum Z = Zb + Z g will be 
dominated by the term with the largest value of q. Thus, the point P has central charge 

/ 6e 2 6e 2 \ 

c(n b , rig) = max 1 - - — ^— , 1 - - — , (33) 
\ 1 - e b 1 - e g J 

and, by the usual identification of the critical Potts model with the dense phase of the 0(n = Jq) model, the critical 
exponents are simply those of a single 0(max(nt>, n g )) model in the dense phase. 

We would expect that only this large- w c portion of the phase diagram gets modified by letting the contact interaction 
be flavor dependent. Let us recall that in the conventional 0(n) model Q with a finite positive vacancy fugacity w c 
the critical behavior of the loops is described by either of two critical branches. The first branch, known as the dense 
branch j2| E8|, is attractive in w c and as such controls the entire domain of low w c . Its central charge is the one 
referred to above: 

c = 1 - 6S 2 /(1 - e) (34) 

in the usual parameterization n = 2cos(-7re). The second branch, known as the dilute branch ^ |27|, is repulsive in 
w c and as such requires w c to be tuned to a particular ri-depcndcnt critical value. In other words, the fugacity of a 
vacancy can tune the O(n) model to its critical point. The central charge of the dilute phase is 

c= l-6£ 2 /(l + e) , (35) 

using the same parameterization as above. 

In particular, in the DPL 2 model |23] the two loop flavors act as decoupled O(n) models, and depending on the 
fugacities of the two flavors of vacancies each of the models can reside in either the dense or the dilute phase, giving 
a total of four different phases. We expect this conclusion to hold true in the generalized eight- vertex model (i.e., 
with an added bending rigidity wx)- Note that only when rib = n g can we simultaneously take the two 0(n) models 
to their critical point by tuning a vacancy fugacity w c which is common for the two loop flavors. In the general 
case, when rib ^ n g we would need two distinct parameters, u> Cj b and w c>g , as indicated on Fig. [| Presumably this 
would lead to a richer phase diagram, with critical lines corresponding to dense-dilute, dilute-dense and dilute-dilute 
behavior of the two O(n) models. 
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Let us return to the phase diagram shown in Fig. ICj. In the special case of the eight-vertex model, the two oblique 
transition lines shown on Fig. [h] are known to be of slope 1/2 fiofl . Actually, they are just images of the line OM 
under certain exact symmetries of the eight- vertex model |10|. Thus, they have again c = 1, whereas the "bulk" of 
the phase diagram is non-critical. 

These two features can be accounted for within the framework of the generalized eight- vertex model with (rib, n g ) = 
(1,1). First, note that our field theory predicts that the region on Fig. [l(] which is limited by the two oblique lines 
and the coordinate axes is actually critical with central charge c = (dense phase); this is obtained by setting 
eb = e g = 1/3 in Eq. (|32|). This is not in contradiction with the exact result flC)|| that this same region is non-critical 
within the eight-vertex model. Namely, the generalized eight-vertex model is embedded in a much larger Hilbert 
space. More precisely, our statement is that the first and third height components possess critical fluctuations and 
constitute a c = theory, even though the second height component is massive. This scenario is brought out very 
clearly by the numerics, as we observe the leading transfer matrix eigenvalues in the sectors determining the free 
energy and xt to coincide within the concerned region. Thus, xt and K22 vanish identically, cf. Eq. (p6|). 

Second, the field theory also accounts for the fact that, within the 8V model, the two oblique lines are critical with 
c = 1. Namely, we claim that they simply correspond to dilute-phase behavior within the generalized 8V model. 
More precisely, since rib = n g the two decoupled O(n) models must be driven to their critical points simultaneously 
by tuning the common parameter w Ci b = w Cyg = w c . Setting e = 1/3 in Eq. ( |35| ) gives a contribution of c = 1/2 for 
each of the models, whence Ctotai = 1/2 + 1/2 = 1 as expected. 

Finally, in the 8V model, the part of the w c axis with < w c < 1 constitutes a further image of the line OM under 
an exact symmetry. We believe this to be "accidental" in the sense that we have seen no sign of a finite interval of 
the w c being critical within the generalized 8V model with other values of the fugacities. 

Taking a common contact parameter u> Ci b = Wc.g for the generalized 8V model with rib 7^ n g destroys the criticality 
of the two oblique lines of Fig. [n]. They still act as transition lines in the sense that they separate the basins of 
attraction of the dense phase and the points P and T respectively. However, the transition is now expected to be 
a first order one. This is confirmed by our numerical results for the (rib,ri g ) = (0,1) case which show that the 
effective central charge develops violent finite-size effects upon approach of the transition lines. Further support for 
this scenario is furnished by Monte Carlo simulations p2| where a finite concentration of empty sites was shown to 
lead to a first order transition. 

The oblique lines in Fig. [l(] are expected to move away from their exactly known 8V positions when we vary the 
loop fugacities away from the trivial values (rib = n g = 1). Some evidence for this is already available from our 
determination of the melting point M in the Flory case; see Eq. fl27|). In general, we have been able to numerically 
determine the position of the uppermost line from the transfer matrix spectra. Recall from the discussion near Eq. ( |2q ) 
that the coupling K22 can be linked to the gap between the leading eigenvalues in two topologically characterized 
transfer matrix sectors. By scanning through w c at fixed wx we have observed (at least in the Flory case) that these 
two eigenvalues become degenerate as soon as w c moves away from zero (even at a value as small as w c ~ 10~ 6 ). This 
degeneracy eventually disappears when there is a level crossing in the groundstate sector of the transfer matrix, at 
some finite w c . 

We have measured the position of this level crossing as a function of system size and extrapolated it to the 
thermodynamic limit. To test the reliability of the method, we have first applied it to the (rib, n g ) = (1, 1) case. Our 
final estimate w c = 1.52 ± 0.02 at fixed wx = 1 is in good agreement with the exact result w c — 3/2 [[l0|. The same 
method applied to the Flory case, (rib, n g ) = (0, 1), yields w c = 1.4294 ± 0.0005 at wx = 1 and w c = 1.958 ± 0.005 at 
wx = 2. These values are clearly different from those predicted by the 8V model. 

In conclusion, we believe that it would be most interesting to study the generalized eight-vertex model in more 
detail, using the exact techniques of integrable systems. In particular, it is conceivable that the present treatment 
misses some subtle exceptional points in the phase diagram. 

VII. DISCUSSION 

The semiflexible loop model was defined as a generalization of the two-flavor fully packed loop model on the square 
lattice, by introducing a vertex weight associated with vertices at which the loop does not undergo a 90° turn. We 
have proposed an effective field theory of the semiflexible loop model based on its height representation. This leads 
to exact results for the Flory model of polymer melting in two dimensions. Furthermore we have shown that the 
loop model provides a generalization of the eight-vertex model with an interesting phase diagram. Here we comment 
further on these two main results. 
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A. Scaling of semifiexible compact polymers 

Polymers configurations are random and as such they are described by probability distributions. Their critical 
nature, in the long chain limit, is revealed by the fact that these distributions have scaling forms characterized by 
universal exponents. The simplest distribution is the probability p{r, I) that the end-to-end distance equals r for a 
polymer of contour length I. In the scaling limit, when r is much greater than the lattice spacing and much less than 
the radius of gyration of the polymer, we have [|l] 

p(r,l)=r e f(r/n ■ (36) 

Here / is a scaling function, v is the " swelling exponent", and 9 the "cyclization exponent". 

For semifiexible compact polymers, which correspond to the Flory model with wx < w x , the swelling exponent 
is v = 1/2. This is an exact result which simply follows from the fact that compact polymers are space filling, 
regardless of wx- Furthermore, the swelling exponent can be related to the string dimension x 2 ,o through the scaling 
law v — (2 — X2,o) _1 [^5|- Then replacing eb = 1/2 and e g = 1/3 in Eq. (|2l]) gives X2,o — and v = 1/2, for all values 
of wx- This calculation then serves as a non-trivial check on the field theory. 

The cyclization exponent is related to the scaling dimension associated with one real and one ghost loop segment: 
9 = —2x1,1 jisfl - From Eq. ( p3| ) it follows that 9 will vary continuously as the polymer is made stiffer by increasing 
wx- At the melting transition the exact result for 9 follows from the computed value of the critical elastic constant, 
Eq. ©, 




The negative value implies that at the transition (and slightly below it) the two ends of the polymer feel an effective 
attraction. This is surprising as the naive expectation is that the two ends of a polymer will feel an effective repulsion 
due to the self- avoiding constraint. For stiff compact polymers this naive expectation is not met. Whether this will 
persist in three dimensions is an interesting open question. 

B. Generalized eight- vertex model 

The generalized 8V model gives a quite detailed modelization of two-dimensional lattice polymers. It possesses the 
following features: 

• Steric constraints (self-avoidance and connectedness of the polymer chains) are modelled exactly; 

• Possibility of introducing polydispersity, by taking rib away from zero; 

• A bending rigidity parameter wx allows to control the transition between a melt and a crystalline phase; 

• A contact interaction parameter w c (or alternatively a fugacity of a vacancy) controls the transitions between 
compact, dense, dilute, Potts-like and non-critical phases; 

• Possibility of introducing non-local interactions (although of a peculiar form), by taking n g away from one. 

The phase diagram of a somewhat similar model was studied in the Bethe approximation by Lise, Maritan, and 
Pelizzola |2!| . However, in the compact limit the results of these authors are equivalent to Flory's mean-field treatment, 
as they do not take into account the non-local features of the polymers. We have here treated the excluded-volume 
effects in an exact manner. On the other hand, the model of Ref. |^] includes an additional feature: 

• A contact interaction between non-consecutive monomers that are nearest neighbors on the lattice allows to 
drive the model to tricriticality, i.e., to access the theta-point physics. 

This interaction is not present in the generalized 8 V- model. If we were to include it, we would need the contact 
interaction to be flavor dependent (the authors of Ref. [^9] do not consider what we refer to as ghost loops). However, 
we do not believe that anything new can be learnt from such a generalization. First, the contact interaction is 
redundant in the compact phase (w c = 0), as the number of contacts is constant (actually maximal) in any fully 
packed configuration. Second, in the non-compact phases (w c ^ 0) our field theory predicts a decoupling into two 
independent O(n) models. One would expect the flavor-dependent contact interactions to act independently on the 
two decoupled models, and the problem essentially reduces to that of the theta-point physics of a standard O(n) 
model @. 
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FIG. 11: Reference state used for generating the sector with flavor-A; strings (k — b,g). 

MUW" 

FIG. 12: Reference state used for generating the thermal sector. 

We leave it as an interesting question whether the generalized 8V model can be tackled using the methods of 
integrable systems. From Fig. |^ it can obviously be formulated as a forty-vertex model (taking into account the loop 
orientations) with complex vertex weights. To our knowledge, such a model has not been studied previously. If one 
could solve it, it would be particularly interesting to work out the exact expression of the coupling constant K-2i{wx) 
as a function of the loop fugacities rib and n g . Also, it is conceivable that our field theoretical approach has missed 
some exceptional critical points in the phase diagram. 



C. Order of the melting transition 



In this paper we have established that the order of the melting transition within the Flory model is second order, 
as first suggested by Saleur 1 13 . We have also explained how the introduction of a finite density of vacancies may 
lead to a first order transition, as observed in Monte Carlo simulations This combined scenario settles a long 

controversy in the literature || . 



APPENDIX A: CONSTRUCTION OF THE TRANSFER MATRICES 

The transfer matrix construction of Ref. jl^] relied on an explicit bijection between the set of allowed connectivity 
states C and the set of integers Z\r\ — {1, 2, 3, ... , \C\}. However, in many cases it is difficult to furnish an a priori 
characterization of the set of allowed basis states and its cardinality. Moreover, some of the states utilized in ju| were 
found to carry zero weight in the leading eigenvectors of the corresponding sectors of the transfer matrix, and so one 
should think that it would be possible to eliminate them from the outset. 

To remedy this situation it is preferable to use another approach. Without prior knowledge of the state space, the 
latter is explicitly generated by acting with the transfer matrix T on a reference state \vo) which is known to belong 
to the image of T in the concerned sector. In this way, a certain number of image states is generated, which can 
be inserted in an appropriate data structure using hashing techniques (3lj. One then acts with T on these states, 
generating a new list of states, and continues in this way until no new states are generated. The resulting list is the 
complete state space of T in the concerned sector. 

It remains to find an appropriate reference state \vq) for each physically interesting sector of T. For the sector 
(sb, s g ) in which Sk flavor-A; strings (k = b, g) span the length of the cylinder generated upon action of T, the reference 
state can be chosen as shown in Fig. [ll]. This state simply consists of (L — Sb — s s )/2 real arches followed by Sb real 
strings and s g ghost strings. L must of course have the same parity as Sb + s s . 

Choosing the sector corresponding to the thermal scaling dimension is a little less obvious. A useful observation is 
made by exploiting the correspondence with states of the six-vertex model, as depicted in Fig. [| In a given row, let 
jyeven &n( ^ ]\[° dd b e the number of flavor-A: loop segments residing on even and odd vertical edges, respectively. Then 
define 

Q = (N° vcn - N° dd ) - (7V g cvcn - 7V g odd ) . (Al) 

By inspection of Fig. || it is seen that Q is nothing but the vertical flux of arrows within a given row. By the ice rule, 
Q is a conserved quantity and can thus be used to label a sector of T. 

The reference state of Fig. [ll] with (sb, %) = (0, 0) is seen to have Q = 0. The first excited state with no strings 
has Q = ±4 and is depicted on Fig. O] Its first four sites are occupied by mutually penetrating real and ghost arches, 
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1720 


196 


684 


1056 


10 


2784 


1280 


6960 


21760 


2520 


10320 


14800 


12 


25888 


13605 


77664 


283584 


32565 


151500 


205920 


14 


259382 


149604 


907837 




425019 







TABLE II: Sizes of various sectors of the SFL model transfer matrix defined on a cylinder of width L with periodic boundary 
conditions. The symbol (sb, s g ) labels the sector in which Sk flavor-fc strings (k — b,g) run along the length of the cylinder. 
The sectors (0, 0) and 'Thermal' have no strings, but the parity of the number of flavor-crossings in the basis states is fixed to 
be even and odd, respectively. 



followed by (L — 4)/2 simple real arches. In general, for any given L, states with Q — ±4g exist for q = 0, 1, . . . , [i/4j . 
The number of states in the k'th sector is just 

L/2-q 



E ( L ( 2 )( L/2 )c L/ ^c n , m 

z — ' \n + qj \n - qj 



n—q 



where C n = ^r^jjrni are the Catalan numbers. Using a sum rule on the binomial coefficients, it is easily seen that the 
total number of states without strings, summed over the sector index q, reads simply 



L/2 , L 



2/i 



E 



CL/2-nCn- (A3) 



This is nothing but the dimension of the state space used in 15 . 

From entropic reasons it is fairly obvious that the free energy belongs to the sector q = 0. We are now going to 
argue that the thermal scaling exponent is linked to the gap between the first eigenvalue in the q = and 5=1 
sectors, cf. Eq. (E51). The first reason is that, by cons truc tion, the q = 1 constraint acts as an excitation within the full 



state space (with all values of q included, as in Eq. (A3)), and hence should correspond to a subdominant eigenvalue 
within that space. Indeed, it is observed numerically that the second eigenvalue obtained from the transfer matrix of 
Ref. jl5) coincides with the leading eigenvalue of the q = 1 sector, obtained by using the techniques outlined above. 

As a second argument, note that in the language of the SFL model height mapping, encircling the first four sites 
of Fig. [l2| yields a height dislocation of A C + B — D. By the four-coloring rule, A-t-B + C + D = 0, this is the 
same as 

m T = 2(A + B) = -2(C + D). (A4) 

But the latter is also the height defect associated with a defect vertex (C, D, C, D) that corresponds to excluding the 
real loops from that vertex, which is exactly a thermal-type excitation (and to wit the one that is used for computing 
the critical exponent xt within the field theory). 

In the field theory, one might compute the exponent corresponding to the insertion of magnetic defects ig'mj- at 
either end of the cylinder. In the transfer matrix, these should simply be linked to the gap between the sectors q = 
and q = q' . 

In Table we show the sizes of the various transfer matrices used in this work. The columns labelled (0, 0) and 



Thermal' correspond to the expressions (A2) with q = and q = 1, respectively. For the other columns, similar 



arc 



expressions may be worked out along the lines of |15[ . 

Finally, let us remark that the computations for the generalized eight-vertex model introduced in Section VI B 
produced from the same reference states, but slightly generalizing the transfer matrix to accommodate the contact-type 
vertices shown on Fig. KJ. 
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We shall usually consider rib and n g real and non-negative, though we believe the model to have interesting properties for 
negative, or even complex, fugacities. 

The vertex weights of the six-vertex model are conventionally called a, b and c. We shall henceforth trade the variable c 
for «jx to avoid confusion with the central charge, which will play a major role in what follows. 

Note in particular that the field theory construction only applies to geometries that are isomorphic to the sphere with 
a point at infinity under a conformal transformation. Generalizations to higher genera have so far only been treated in 
mean-field theory |(J. 



